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The circular and Jacobi ensembles of random matrices have their eigenvalue support on the 
unit circle of the complex plane and the interval (0, 1) of the real line respectively. The 
^ • . averaged value of the modulus of the corresponding characteristic polynomial raised to the 

I power 2/i diverges, for 2/z < —1, at points approaching the eigenvalue support. Using the 

I theory of generalized hypergeometric functions based on Jack polynomials, the functional 

form of the leading asymptotic behaviour is established rigorously. In the circular ensemble 
■ case this confirms a conjecture of Berry and Keating. 



m : 1 Introduction 

Random matrices from the classical groups — N x N unitary matrices U(N), N x N real orthogonal 
matrices 0{N), and N x N unitary matrices with real quaternion elements embedded as 2N x 2N 
' complex unitary matrices Sp{N) — play a special role in the application of random matrix theory to 

. number theory (see e.g. the recent review Of particular interest in such applications are the random 

, matrix averages of the modulus of the characteristic polynomial raised to some power 2/i say. Thus in 

! HO] it was shown how knowledge of 

(|det(z/-i7)|2^\ (l.I) 
CLi' \' ^ ^' /ueu{N) ^ ' 

' for \z\ = 1 (in this case Hl.l|l is in fact independent of z) allows the mean value of the 2/i-th power of the 

' modulus of the Riemann zeta function on the critical line to be predicted. Knowledge of the analogue 

of for the classical groups 0{N) and Sp{N) allows for similar predictions in the case of families of 
L-functions |S1II3E1- 

\ Our interest is in the asymptotic behaviour of Ijl.lf) and its analogues as\z\ approaches unity. Consider 

' in particular the generalization of Hl.l|l 

N 

TTlz-e^^f^) , (1.2) 

1=1 

where C/3Ejv (circular /3-ensemble) refers to the eigenvalue probability density function proportional to 

Yl |e^9. _g»e,|/3^ -Tr<ei<Tr. (1.3) 

l<j<k<N 

The cases (3—1 and P — 4 are well known in random matrix theory as the COE (circular orthogonal 
ensemble) and CSE (circular symplectic ensemble) respectively, while the case /3 = 2 is just (|1.1|) . For 
2/i < —1 it has been argued by Berry and Keating fSI that 1)1. 2|l diverges as |1 — \z\\^^ where, provided 
^ |/i|, (5 as a function of /i is the non-smooth function 

S = int[^^ + 1] (fc - 1 + ^ - f + 1]) ' ^ ^ 21^1' (^-^^ 

with int[-] denoting the integer part. It was also suggested that for k values such that (fc— 1) //3 is an integer, 
there is a logarithmic correction to the behaviour of 1)1.2(1 which then diverges as |1 — \z\\^^ log |1 — |z||. 
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The intricate behaviour exhibited by H1.4|l is due to a phenomenom termed singularity dominated strong 
fluctuations. Thus for 2/i < — 1 and \z\ close to unity, degeneracies of the spectrum of varying order are 
being probed A similar mechanism has been identified in the analysis of twinkling starlight PP, van 
Hove- type singularities and the influence of classical periodic orbit bifurcations on quantum energy 
level 0] and wavefunction statistics |18j . 

The arguments presented in [3] leading to (|1.4|) were heuristic. Subsequently Fyodorov and Keating 
[TT] gave a rigorous derivation for the exponent (|1.4|l in the case /3 = 1, 2/i e Z<o, for the Gaussian 
orthogonal ensemble analogue of (|1.2|) , 



det(a;/-X)P'^\ . (1.5) 

/XGGOEjs 



Thus they proved that for N ^ ii and k = 2|/i| e Z>o 

\ / XgGOEjv e-+o e 

(the explicit evaluation of Cnm was also given). In this paper we will give a rigorous derivation of the 
exponent (|1.4(l . together with the logarithmic corrections when present, for the average (|1.2(l in the case 
of general f3 > and general 2/z < —1. Our analysis relies on identifying (|1.2|l as a special generalized 
hypergeometric function based on Jack polynomials j25[ I17| . We are then able to use known asymptotic 
properties of the latter to deduce the sought asymptotic behaviour of H1.2|l . 
Also studied will be the x ^ 1+ asymptotic behaviour of 

N 

(ni--H>)„,. (1.7) 

where J/JE^v (Jacobi /3-ensemble) refers to the eigenvalue probability density function proportional to 

N 

Ylxlil^xif II \xk-x,f, 0<xi<l. (1.8) 

1=1 l<j<k<N 

In the special case /3 = 2, (|1.7(l includes the analogue of (jLlll in the cases of the classical groups 0{2N+1), 
0{2N) and Sp{N). Again our analysis proceeds via identifying H1.7|l as a generalized hypergeometric 
function based on Jack polynomials. 

As a final issue, it is pointed out in ^21 that for fixed \z\ ^ 1 and N ^ oo, Szego's theorem |23 on 
the asymptotic behaviour of Toeplitz determinants with smooth generating functions implies (|l.l|l is a 
simple Gaussian in /i, 

lim /|dct(z/-[/)|2f\ ^ g-A.= iog|i-|.|=|^ 1^1^^^ /^ 9) 

N^oo \ I U£U{N) 

This gives for the exponent characterizing the divergence as \z\ — > 1~, 6 — and so does not reproduce 
(|1.4|l . Therefore the limit N ^ oo with |z| 7^ 1 fixed has the effect of smoothing out the singularity 
dominated strong fiuctuations. Note that in this limit the point z is a macroscopic distance away from 
the eigenvalue support, as measured in units of the inter-eigenvalue spacing. A theorem of Johansson 
^lEI will be used to establish the analogue of (|1.9(l for the average (ll.2|l . 
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2 Generalized hypergeometric functions based on Jack polyno- 
mials 



2.1 Definitions 

The Gauss hypergeometric function 2F1 is, for < 1, defined by the series 

2Fi{a, h- c; x) ^ ^^j^^^ ■ (2.1) 

A muhivariable generahzation of (|2.1I) can be defined in which the sum over the non-negative integers k 
is replaced by a sum over partitions k :— . . . , kat), > Kj {i < j) and G Z>o; the function 

iu)k:^uiu + l).--iu + k-l)^^^^-^ 

T{u) 

is replaced by 



N 

r{u - (j - l)/a) ' 



fc! is replaced by \k\\ :— (X^jLi'^j)-! ^^'^ replaced by the homogeneous symmetric polynomial 

ci°'\xi, . . . ,xn) which in turn is proportional to the Jack polynomial Pk°'\xi, . . . ,xjv)- Thus for |a;i| < 
1, . . . , I^Arj < 1 one defines PSIITTI 

(a, 6; c; XI, ... , x^) := ^ M|^ci") (xi , . . . , (2.2) 

Let us say a little on the definition of the Jack polynomials Pk°'\xi, . . . ,xn) =■ Pk°'\x) and their 
renormalized version ci"\xi,...,XN) =■■ Ci'"\x). The former are the unique homogeneous polynomials 
of degree |k| with the structure 

Pi")(z) = m« + ^a,^m^ (2.3) 
(the a^fj, are some coefficients in Q(q;)) and which satisfy the orthogonality 

where 

/•1/2 /•1/2 

•^-1/2 J-1/2 l<j<k<N 

In (|2.3I) denotes the monomial symmetric function (polynomial) in xi, . . . , arjv corresponding to k, and 
II < K means '^j — X]j=i each p — 1, . . . , TV. Regarding the definition of ci"^ {xi , . . . , xjv), 

let 

rfK= n («(a(«,j) + l) + ^(*,j)), (2.4) 

where the notation (i, j) € k refers to the diagram of k, in which each part Ki becomes the nodes («, j), 
1 ^ J ^ '^i on a square lattice labelled as is conventional for a matrix. The quantity a(i,j) is the so 
called arm length (the number of nodes in row i to the right of column j), while is the leg length 

(number of nodes in column j below row i). In terms of d'^ we have 

C^^\z):^^^P^-Hz). (2.5) 
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2.2 Integration formulas 

Introduce the Selberg integral 



■^0 Jo ,_n n^,-,,.^Ar 



1=1 l<j<k<N 



N 



S'jv(Ai, A2, A) 
and the Morris integral 

.1/2 .i/Z 

■^^1/2 "'-1/2 J^;^ l<j<fc<JV 

(both these integrals have gamma function evaluations — see e.g. [7| — although these will not be needed 
here) . Fundamental to our ability to express p.7|l and (|1.2|l in terms of the generalized hypergeometric 
function ()2.2|) are the generalized Selberg integral evaluation |52 ^1 [T7] 

N 



Sn{X 



1 .1 .1 

l,A2,l/ajJo Jo l^.tt^^r 

(2.6) 



pM^iv^ [Ai + (iV-l)/a + l] 

[Ai+A2 + 2(7V-l)/a + 2]l") 



1 l<j<k<N 
(a) 



.1/2 .1/2 W 

J-1/2 J-1/2 fj^^ 
X [] |g2..e. _g2.»e,|2/« ^p^a)(^W) t3!ll (2.7) 



and the generalized Morris integral evaluation [Hj 

»l/2 .1/2 N 

do • ■ • / 

Mjv(a,fo,l/a) 7-1/2 ^ 7-1/2 

i<,<;c<A' [l + a+(iV-l)/a](") 

f H2.6|l and H2.7|l are in fact equivalent JS]), where 

a;i— a^iv — 1 

Equally important is the generalized binomial summation formula [251 123| 

r ^ 

,Ft\a-, x,,...,x^):=J2 TTT^^i^H^i. . . . ,x^) ^ - x,)-\ (2.8) 

Thus we see by multiplying both sides of H2.6|l and H2.7|l by 

(xa)l"l[r]i") 
d'. 

and summing over k, using H2.8(l on the left hand sides and using the definition H2.2|l on the right hand 
sides together with the fact that is homogeneous of degree |k|, that |H] 
^ .1 .1 Af 



S'Ar(Ai, A2, l/a) Jq 



.1 .1 " 



1=1 l<j<k<N 

1 2 I 

,X'"\r, -{N - 1) + Ai + 1; -(iV - 1) + Ai + A2 + 2;ti, . . . , i^) (2.9) 
a a \ti=---=tN=x 



and 



nl/2 .1/2 N 

— -— T — , . / dOi ■ ■ ■ t 

MN{a,b,l/a) J_i/2 J ^1/2 

n -e2"«^p/" = 2^i"^(^,-&;-(^-l)+a + l;ti,...,ijv) (?.10) 



J-1/2 J-1/2 , , 



X 

l<j<fc<Ar 
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3 Random matrix averages as generalized hypergeometric func- 
tions 

Here we will express (ll.2|l and (|1.7|l in terms of 2-^1"'' • Consider first (|1.2|) . By definition and simple 
manipulation 



N 



5' ■ 

1/2 /.I/2 N 

d0N 

1 = 1 l<j<k<N 



Now 

Regarding the integrals over 9i G [—1/2, 1/2] as parametrizing the unit circles wi — e^'^'^' in the complex 
plane, then deforming each wi by wi 1-^ ze^'^*^' (this is immediate for 2/i, /3/2 G Z>o by Cauchy's theorem; 
it remains valid for general complex /i and /3 by Carlson's theorem) we see that we can write 

N 



n 

1=1 



Mn{0 



C/3E„ 

P/2j J_i/2 J-1/2 77 



'-1/2 J-l/^ i=i 

l<j</£<A' 

f-1/2 ^1/2 JV 



MAr(0 



■i— 7— / d9i--- f d6'Ar TT|l + |zpe2"«'P'e-"^'''|l + e-2"Si|M 
, 0, p/2)7_i/2 J-1/2 



X II |e2^»e, _g2«e,|/3 



n 

l<j</£<Ar 
,(2//3)^_,, 

2 

where the third equality follows from (|2.1U|) . and is valid for |zp < 1. For |zp > 1 one notes the simple 
symmetry 

N N 



2Fr'^^(-M,-M;^(A^-l) + l;ii,...,t^) , ,^ (3.1) 



e ' 

C/3Ejv ■ ■ '^f'J^'^ /C/3Ej, 



thus relating this case back to the case |zp < 1. We remark that in the case P ~ 2 the same generalized 
hypergeometric function evaluation H3.1|l has recently been given in |10| . 

Consider next (|1.7(l . By definition and simple manipulation, for x real and greater than or equal to 
unity 



AT 



J/3E 



/ ■ • • / dtNW ttii - ti)\i - ti/xf^ n \tk - 1 



SNia,b,(3/2) 







1 = 1 l<j<k<N 



^2^N 2j^^(2//3)(_2/i, Un -l) + a-l3{N -l) + a + b + 2;h,...,tN)\ (3.2) 
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where the second equahty follows from (|2.9|l . The task now is to analyze the generalized hypergeometric 
functions in and H3.2I) as \z\ 1 and x 1+ respectively. Fortunately, as will be presented in the 
next section, the required asymptotic formulas are available in the literature. 

4 Asymptotic forms 

Using the series definition (|2.2|) . Yan has analyzed the behaviour of the function 2 -F'l"'' asxi, . . . ,xjv 
1^ . Thus according to Proposition 4.4 of |2S1 the following result holds, where the notation A{x) « B{x) 
means that there exists two positive numbers Ci and C'2 independent of x such that 

B{x) 



Proposition 1. Let 

and suppose for all k 



We have for ~1 < Xi < I (i = 1, . . . , A^) 
(i) If"f> {N -l)P/2, then 



J := a + b — c. 



> 0. 



N 



2F^"\a, b;c;xi,...,XN)- Y[{1 - x,)-\ (4.3) 



(a) If ^ < —{N — 1)13/2, then there exists a constant C such that 

2F^''\a,b;c;xi,...,XN)<C. (4.4) 
(Hi) Ifj = + j - 1), j = 1, . . . ,iV, then for xi = ■ ■ ■ = xn = t, 

,Fi''\a,b;c;x^,...,XN) « (1 _ t)-0-ibW2 log -i- (4.5) 

Xi — --- — XN—t 1 ~ t 

(iv) IfPi-^ + J - 1) < 7 < /3(-^ + j), i = 1, . . . , iV ™ 1, then for xi ^ ■ ■ ■ ^ xn ^ t, 



,F[''\a,b;c;x„...,x^) « (1 _ i)-.(7+(Jv-,)W2)_ (4g) 

Xi — ---—XN—i 

Application of PropositionQlto the 2-^1"'' evaluations H3.1|l and H3.2() gives the desired functional form 
of the singular behaviour of p.2|l and p.7|l in the case of a negative exponent 2/i. 

Corollary 1. Consider the circular ensemble average jJ.g)) . Suppose /i < and |z| < 1. For 

2|m| = /3(j - 1) + 1 {j = l,...,N) 

we have 

N 

(ni-»*l"'l,, -a-W)-'-«"log^ (4.7) 



1=1 
while for 

/3(j-l) + l<2|^| <^j + l (j = l,...,iV) 
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and thus 

J = int[(2lAi| - + 1], (2|/i| - l)//3 + 1 ^ Z 

AT 

(TT|z-e*«f^) «(i_|^|)-K2|H-i-0-iW2)_ (4 8) 

\ / C/3Ejv 

Consider the Jacobi ensemble average jj. ?[ ). Suppose fi < and a; > 1. For 

= /3(j - 1) + 2 + 6 



have 



while for 



and thus 



we have 



N 

-^^^(H \x ~ xM « (1 - log — (4.9) 



/3(j-l) + 6 + 2<2|Ai| </?j + 6 + 2 (j-l,...,iV) 
j = int[(2|/i| - b - 2) IP + 1], (2|/i| - & - 2)//3 + 1 ^ Z 



a; 



-2/j7v/ 



J/3E 



The asymptotic behaviour H4.8|l indeed exhibits the exponent l|1.4() predicted by Berry and Keating 
0, and the existence of a logarithmic correction in the cases that (2|/i| — 1)/ (3 is an integer is confirmed. 
Finally, let us apply the results H4.9(l and H4.1()|l to deduce the asymptotic z ^ 1+ behaviours of 



/|det(z/-C/)p'^\ , ^l<Q 

\ I ueG 



forGthe classical groups S'p(iV), C'+(2iV) andO-(2iV) (we leave the cases of 0+(2iV-l) andO-(2Ar-l) 
as an exercise for the interested reader). Consider first Sp{N). The eigenvalues come in complex conjugate 
pairs e^*^^ {j ~ 1, . . . , N) so we have 

N 2,-1 

det(z/ -U) = i2zf Y[ (^-^ - cos 9^ . (4.11) 



Furthermore, we know (see e.g. T) that in the variable cos9j =: Xj the eigenvalue probability density 
function is of the form (|1.8(l with a = b = 1/2, /3 = 2. Thus setting z=:l + e, 0<e^l and noting 
(z^ + l)/2z - 1 + 6^/2 it follows from (O that for 



2|/i|=2j + l/2 (j = l,...,iV) 

we have 

/ldetffl + e)/-C/)P'^^, 

ueSp(N) 



(\detiil + e)I-U)\"^) «6-2(^-iWlogi (4.12) 



while for 

2j + 1/2 < 2\^l\ < 2(j + 1) + 1/2 

we have 

/|det((l + e)/-C/)p^\ ^g-2j(2|/.h3/2-j)_ (4^^3) 
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For the classical group 0^{2N) the eigenvalues again come in complex conjugate pairs e^*^^ [j — 
1, . . . , iV), so (|4.11|l remains valid. In the variable cos6'j —: Xj the eigenvalue probability density function 
is proportional to (|1.8|l with a = b = -1/2, /3 = 2 (see e.g. |7])- Thus setting z = l + e, 0<e<lif 
follows from H4.9|) that for 

2|m|=2j-1/2 (j = l,...,iV) 



we have 

/|det((l + e)/-f/)pf\ «e-20-ibiogi (4.14) 

\ /ueo+(2N) e 



while for 

2j-l/2<2|^|<2(j + l)-l/2 

we have 

/|det((H-e)/-i7)|2^\ ^g-2j(2|M|-i/2-j)^ (4^^5) 

\ /ueO+(2N) 

In the case of the classical group 0~ {^N), there is a pair of fixed eigenvalues at ±1, with the remaining 
eigenvalues coming in complex conjugate pairs e**^^ (j = 1, . . . , A*" — 1). Hence in this case 

^'-i 2 I -1 

det(z/-C/) = (^2_i)(2z)^-i Yl (i-^-cos0,). (4.16) 



In the variable cos^j =: xj the eigenvalue probability density function is proportional to H1.8(l with 
N ^ N - I and a = 1/2, 6 = -1/2, /3 = 2. As only the value of b enters in and ^J^, the sole 

modification of H4.14|l and (|4.15() is multiplication by e to account for the factor of — 1 in (|4.16() . Thus 
for 

2|/i| = 2j-l/2 (j-l,...,iV-l) 



we have 



while for 



we have 



/|det((l + e)/-C/)|2'^\ « e-20-i)J+i log - (4.17) 
\ /ueo-{2N) e 

2j-l/2<2|^| <2(j + l)-l/2 (j = l,...,7V-l) 

/|det((l + e)/-;7)pA'\ ^g-2,(2|H-i/2-j)+i. (4.18) 
\ /ueo-(2N) 



5 The macroscopic limit 

Consider the general Toeplitz determinant 



r 1 /•27r 

I5Ar[e"(^)] := det — / e^^^^e'^^-'^)^ ^6* . (5.1) 



2^ Jo 



j,k=l,...,N 



Such determinants are related to averages over U{N) by the simple to establish formula 

N 

DnIc^^'^] = (Y\ e"^''^) ■ (5-2) 

\ / U(N) 
1=1 ^ ' 

With a{9)) = J2^=-oo flp6*P^ and the Fourier coefficients falling off fast enough that 

00 

J2 \P\K\^<^^ (5-3) 
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an asymptotic formula of Szego gives 

oo 

DNle"^'^^] - expf7Vao + y"papa-p + o(l)y (5.4) 

p=l 

Now, in the case 

a(6l) = log |z-e''^ I, \z\ < I, (5.5) 
the Fourier coefficients have the exphcit form 

ap = -^(p>0), ap = ^^(p<0), an = 0. (5.6) 

Since we are assuming < 1 these coefficients satisfy (|5.3f) . Substituting H5.6|) in (|5.4ll gives (|1.9II . 

In this section, as our final issue, we will make use of a generalization of (|5.4I) due to Johansson to 
establish the generalization of (|1.9() when the average over U (N) is replaced by an average over C/SEtv- But 
before doing so, following we make note of the interpretation of p.9|l as specifying the distribution 
of the linear statistic 

N 

A(z)=^log|z-e^^^■|^ |z|<l. (5.7) 

Thus let Pz{t) denote the probability density that A(z) takes on the value t after averaging over the 
eigenvalue distribution of U{N). Then we have 

P.{t) = {s{t-A{z)))^^^^ (5.8) 

and consequently 

Pz(k):^ P,(i)e**'^di = ([|e2''='°sk-e'«M\ . (5.9) 

Thus Pz{k) is precisely (15. 2|) with the substitution (j5.5|l . and fi therein set equal to ik. Consequently 

hm P,(t)=e-'='-'/2^ a2 = -log|l-|zp|2, (5.10) 

telling us that the limiting distribution of (|5.7I) is a Gaussian with an 0(1) variance taking the explicit 
value - log |1 — \z\^\'^. 

To study l|5.8|l with the average over U{N) replaced by the probability density function C/SEn, we 
make note of the following generalization of the Szego theorem (15 .411 due to Johansson EI ■ 

Proposition 2. Suppose a{9) — '}2,^=-oo ^pS*^^ '^^'^ Fourier coefficients Jail off fast enough that the 
inequality i|5.5]) holds. We have 

N oo 

TTe"^''^) - exp fiVao + -5 Vpapa_p + 0(1)). (5.11) 

l=\ p=l 

As an immediate corollary, by substituting H5.5|l for the Fourier coefficients the sought generalization of 
(|1.9|l can be deduced. 

Corollary 2. For \z\ < 1 we have 

N 

lim /TT|z-e^«f^\ ^g-(2MV/3)iog|i-kPl. (5.12) 

N-,oo \f-J^ I C/3En 
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Of course the result H5.12|l affords an interpretation as a fluctuation formula as in H5.10|l . the only 
difference being that cr^ — — (2//?)log|l - 

An obvious question is to seek the analogue of (|5.12|l for the average (|1.7|l with x ^ [0, 1]. Johansson 
[n| has derived the analogue of H5.11|l . but only in the case (3 — 2 and further with the restriction that 
a{9) is a polynomial in 6', which is not the case for H5.5|l . 
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